Solvation and Dissociation in Weakly Ionized Polyelectrolytes 
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We present a Ginzburg-Landau theory of inhomogeneous polyelectrolytes with a polar solvent. 
First, we take into account the molecular (solvation) interaction among the ions, the charged 
monomers, the uncharged monomers, and the solvent molecules, together with the electrostatic 
interaction with a composition-dependent dielectric constant. Second, we treat the degree of ioniza- 
tion as a fluctuating variable dependent on the local electric potential. With these two ingredients 
included, our results are as follows, (i) We derive a mass reaction law and a general expression 
for the surface tension, (ii) We calculate the structure factor of the composition fluctuations as a 
function of various parameters of the molecular interactions, which provides a general criterion of 
the formation of mesophases. (iii) We numerically examine some typical examples of interfaces and 
mesophase structures, which strongly depend on the molecular interaction parameters. 



I. INTRODUCTION 

Polyelectrolytes are much more complex than low- 
molecular-weight electrolytes and neutral polymers^ 2 -! 3 -. 
Above all, the electrostatic interaction among the ion- 
ized monomers on the polymer chains and the mobile 
ions strongly influence the chain conformations and the 
mesophase formation 2 .! 3 .!!. Second, the dissociation (or 
ionization) on the chains should be treated as a chem- 
ical reaction in many polyelectrolytes containing weak 
acidic monomers^! 7 -, which is under the influence of the 
local electric potential. Then the degree of ionization is 
a space-dependent annealed variable in inhomogeneous 
polyelectrolytes, while it has mostly been treated to be a 
given constant in the theoretical literature. Such ioniza- 
tion heterogeneity should not be negligible in structure 
formations and in phase separation. Third, for water-like 
solvents with large dielectric constant e, polymers are of- 
ten hydrophobic and small ions are hydrophilic 6 , which 
can also affect the phase behavior. However, not enough 
attention has yet been paid on the effects of such short- 
range molecular interactions, where particularly relevant 
is the solvation (ion-dipole) interaction between ions and 
polar molecules^. 

In this paper we hence treat the degree of ionization as 
a fluctuating variable and include the molecular interac- 
tions. We show their relevance in polyelectrolytes in the 
simplest theoretical scheme. That is, we use the so-called 
random phase approximation 1 in the Flory-Huggins 
scheme of weakly charged polyelectrolytes 9 - ' 10 ' 11 ' 12 . On 
the basis of a recent Ginzburg-Landau theory of ion dis- 
tributions in binary mixture o 13 ' 14 ' 15 ' 16 , we account for 
the solvation (hydration in aqueous solutions) between 
the charged particles (ions and ionized monomers) and 
the solvent molecules, whose free energy contribution 
usually much exceeds the thermal energy T (per charged 
particle) 8 . Hereafter we set the Boltzmann constant 
equal to unity. 

In one-phase states of weakly charged polyelectrolytes, 
the structure factor S(q) of the composition fluctuations 
with wave number q was calculated in the random phase 



neglected. Let <f> and / be the polymer volume fraction 
and the fraction of charged monomers on the chains, re- 
spectively. In this approximation the inverse of S(q) is 
expressed as 



1 =m+ 



S(q) 



\2aMl - 



Vq(<J 2 + K 2 ) ' 



(1.1) 



approximation 
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where the solvation interaction was 



where the first term f depends on <p and T, vq = a is the 
volume of a monomer, Is = e 2 /eT is the Bjerrum length, 
and k is the Debye-Hiickel wave number. Due to the last 
electrostatic term in Eq.(l.l), S(q) can have a peak at 
an intermediate wave number q m for small k or for low 
salt concentration. We here mention some related exper- 
iments, (i) Such a peak has been observed by scatter- 
ing in one phase states of charged polymer system s 17 ' 18 . 
It indicates formation of mesophases in sufficiently poor 
solvent, as was confirmed for a semidilue polyelectrolyte 
solution 19 , (ii) On the other hand, for neutral, polar 
binary mixture s 13 ' 14 (or a mixture of neutral polymers 
and a polar solvent^) with salt near the critical point, 
we calculated S(q) in the same form as in Eq.(l.l). In 
electrolytes, the solvation interaction can strongly affect 
the composition fluctuations particularly near the criti- 
cal point. In fact, in a recent scattering experiment 2 ^, a 
peak at an intermediate wave number has been observed 
in S(q) in a near-critical binary mixture with salt, (iii) 
We also mention a finding of a broad peak in S(q) in 
semi-dilute solutions of neutral polymers in a polar sol- 
vent with sal t 21 ' 22 , where the solvation effect should be 
crucial. Thus we should calculate S(q) in weakly charged 
polyelectrolytes including the solvation interaction. 

We also mention calculations of the interface profiles 
in weakly charged polyelectrolytes in a poor solvent us- 
ing the self-consistent field theory 2 ^!. In these papers, 
however, the solvation interaction was neglected. In poly- 
electrolytes, the solvation interaction should decisively 
determine the charge distributions around an interface, 
as in electrolyte a 14 ' 16 . In addition, the degree of ioniza- 
tion should significantly vary across an interface in poly- 
electrolytes, because the dissociation process strongly de- 
pends on the local counterion density. 
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The organization of this paper is as follows. In Sec- 
tion 2, we will present a Ginzburg-Landau approach ac- 
counting for the molecular interactions and the disso- 
ciation process. We will introduce the grand potential 
and present a theoretical expression for the surface ten- 
sion. In Section 3, we will calculate the composition 
structure factor generalizing the results in the previous 
theories ^ 10 ' 13 ' 14 . In Section 4, we will numerically ex- 
amine the ion distributions around interfaces and in a 
periodic state. 

II. THEORETICAL BACKGROUND 

We suppose weakly charged polymers in a theta or 
poor solvent consisting of a one-component polar fluid. 
We assume / -C 1 to ensure flexibility of the chains. As 
suggested by Borue and Erukimovich^ and by Joanny 
and Leibler— , the random phase approximation can be 
used in concentrated solutions with 

^>fl\t B /afl\ (2.1) 

We consider the semidilute case 4> > iV -1 / 2 , where N is 
the polymerization index. Then the polymers consist of 
blobs with monomer number g = <$T 2 and length £j 
fin ' in the scaling theory^ and the electrostatic energy 
within a blob is estimated as 

e b = T(fg) 2 £ B /£ b = Tfcf>-*l B a-\ (2.2) 

Here the salt density is assumed not to exceed the density 
of the charged monomers n p . Obviously, the condition 
£b < T yields Eq.(2.1). In our case the Debye-Hiickcl 
wave number k oc (ripls) 1 ^ 2 is sufficiently small such 
that n p ^> k 3 holds and a free energy contribution due 
to the charge density fluctuations oc k 3 is negligibl e 25 ' 26 . 

A. Ginzburg-Landau free energy 

For weakly ionized polyelectrolytes, we set up the 
free energy accounting for the molecular interactions 
and the ionization equilibrium. We neglect the image 
interaction ^ 14 ' 27 and the formation of dipole pairs and 
ion cluster a 12 ' 28 . The former is important across an in- 
terface when the dielectric constants of the two phases 
are distinctly different in the dilute limit of the ion den- 
sities (see comments in the last section), while the latter 
comes into play at not small ion densities. 

The volume fractions of the polymer and the solvent 
are written as <j>(r) and 1 — 4>(r), respectively. For sim- 
plicity, we neglect the volume fractions of the ions and as- 
sume that the monomers and the solvent molecules have 
a common volume v$ = a 3 . Then tf> is also the molar com- 
position. The counterion density is written as n c (r). We 
may add salt with cation and anion densities n\(r) and 
n 2 (r), respectively. The ion charges are Zi with i = c, 1, 
and 2. In the monovalent case, for example, we have 



Z c = 1, Z\ = 1, and Z 2 = — 1, respectively. In the con- 
tinuum limit these variables are smooth coarse-grained 
ones on the microscopic level. 

The number of the ionizable monomers on a chain is 
VmN with Vm < 1. In this work the degree of ionization 
£(r) in the range < ( < 1 depends on the surrounding 
conditions and is inhomogeneous. Then the fraction of 
ionized monomers is 

/ - "mC, (2.3) 
and the number density of ionized monomers is 

n P = v o 1 f ( t ) = u o" 1j/ mC^ ( 2 - 4 ) 

where C and <t> are space-dependent. The condition of 
weak ionization / <C 1 is always satisfied for vm *C 1, but 
we need to require £ <C 1 for vm ~ 1. Furthermore, we 
assume that the charge of each ionized group is negative 
and monovalent (or equal to — e) so that the total charge 
density is written as 

p = e ZiJii — en p . (2-5) 

i=c,l,2 

The overall charge neutrality condition is J drp = 0. In 
more detail in the present of salt, we have 

J dr[Z c n c - n p ] = J drlZxni + Z 2 n 2 ] = 0. (2.6) 

The free energy F of our system is the space integral 
of a free energy density /t in the fluid container. We 
assume that fa is of the form, 

fa = / o (0,T) + ^C(0)|V0| 2 + 4^ jE;2 

+T ^ ni[\n(niV ) - 1 + gi<j>] + T(A + g p <j))n p 
i=c,l,2 

+Tuo W[ClnC+ (1 - C)ln(l - C)]- (2.7) 

In the first line, the first term /o is the chemical part in 
the Flory-Huggins form 1 ^, 

fo = -[-£ In <t> + (1 - 4>) ln(l -4>) + X 4>{1 - 4>)], (2.8) 

where x is the interaction parameter dependent on the 
temperature T and its critical value is \ c — (1 + 
N^ 1 / 2 ) 2 /2. in the absence of ions. The second term 
is the gradient part with the composition-dependent 
coefficient 1 -, 

C{<jj) = l/12a^(l- 0), (2.9) 

1/3 

where a = v Q . The third term is the electrostatic free 
energy, where E = — V$ is the electric field. The elec- 
trostatic potential $ satisfies the Poisson equation, 

V ■ e(0)V$ = -47175. (2.10) 
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The dielectric constant e((f>) changes from the solvent 
value £o to the polymer value e v with increasing <j). For 
simplicity, we assume the linear form, 



e(4>) = Eo + e^ 



(2.11) 



hydrophilic ions and gi < for hydrophobic ions. The 
Born thcory^i and the data of the Gibbs transfer free 
energ y 32 ' 33 both suggest that \gi\ mostly much exceeds 
unity and is even larger for multivalent ions such as Ca 2+ . 



where z\ = e p — Eq. In some binary fluid mixtures, the 
linear form of e(4>) has been measured 2 ^. For charged 
gels, Kramarenko et al. pointed out relevance of strong 
composition-dependence of e((f>) in first-order swelling 
transition 3 ^. In the second line of Eq.(2.7), the coupling 
terms Tg^ni (i = c, 1, 2,p) arise from the molecular in- 
teractions among the charged and uncharged particles, 
while TAo is the dissociation energy in the dilute limit 
0^0. In the third line of Eq.(2.7), we give the en- 
tropic free energy of dissociation^^, where Vq 1 vm4> is 
the density of the ionizable monomers. 

Solvation interaction. We have introduced the molec- 
ular interaction terms (oc <^), which will be simply called 
the solvation interaction terms. In low-molecular- weight 
binary mixtures with sal t 13 ' 14 , such terms arise from the 
composition-dependence of the solvation chemical poten- 
tial of ions pf l , where i represents the ion species. The 
original Born theory 3 ^ gave pf* 1 = Zfe 2 /2e(<fi)Ri, where 
Zie is the ion charge and Ri is called the Born radius. 
Here Ri is of order lA for small metallic ions in aqueous 
solution^. For each ion species, the difference of pf in 
the coexisting two phases is the Gibbs transfer free energy 
typically much larger than T per ion (~ 15T for monova- 
lent ions) in electrochemistr y 32 ' 33 . In polymer solutions, 
the origin of these terms can be more comple x 21 ' 22 ' 34 . 
For example, ions interact with the dipoles of the sol- 
vent molecules and those on the chains differently, affect- 
ing the hydrogen bonding around the chains. Therefore, 
the solvation chemical potential of an charged particle 
of the species i (i = c, 1, 2,p) arises from the interaction 
with the solvent molecules and that with the uncharged 
monomers as 



e|(l- 



(2.12) 



where e™ 1 and ef are the interaction energies. The solva- 
tion contribution to the free energy is given by the space 
integral of the sum, 



Esol 
Pi fli 



= T^2g i (/)n l 



(2.13) 



Here we find gi appearing in Eq.(2.7) expressed as 
9, = Of - "t)/T. 



(2.14) 



The last term on the right hand side of Eq.(2.12) con- 
tributes to a constant (irrelevant) chemical potential for 
i = 1,2 and to the constant Ao in Eq.(2.7) for i = c,p 
(with the aid of Eq.(2.6) for i — c). As a result, we have 
attraction for g^ < and repulsion for g± > between 
the ions (i = c, 1,2) and the polymer chains, while we 
have a composition-dependent dissociation constant from 
i — c,p (see Eq.(2.27)). For water-like solvent, gi > for 



B. Equilibrium relations 

As a typical experimental geometry, our fluid system 
is inserted between two parallel metal plates with area S 
and separation distance L(much shorter than the lateral 
dimension S 1 / 2 ). If the surface charge densities at the 
upper and lower plates are fixed at ±<to, the electrostatic 
energy F e = J dre(4>)E 2 /8n is a functional of 4> and p- 
The potential values at the two plates are laterally ho- 
mogeneous, but are fluctuating quantities 1 ^. 

For small variations S<fi and Sp superimposed on <f> and 
p, the incremental change of F e is written asi 4 - 



SF e 



dr 



e x E z 
8tt 



(2.15) 



Under the charge neutrality J drp — 0, we first minimize 
F with respect to £ (or n p ) and n c at fixed m, n2, and <j>. 
It is convenient to introduce G = F+(TX/e) J drp, where 
TA/e is the Lagrange multiplier independent of space. 
Using Eq.(2.15) we may calculate 8G/8C, and 5G/5n c . 
Setting them equal to zero, we obtain 



c 



= exp[l7 + A-Ao-g p $, (2.16) 

(2.17) 



i-C 

n c = Vq 1 exp[-U - A - g c <f>], 



where we introduce the normalized electric potential by 



U = e$/T. 



(2.18) 



Next, homogeneity of the ion chemical potentials pi = 
SF/Sn, (i = 1,2) yields 



= r#exp[-gi<l>-Zi(U + \)], 



(2.19) 



where nf = u _1 e Mi / T are constants. Here U and A appear 
in the combination U + A in all the physical quantities. 
If U + A is redefined as U, A may be set equal to zero 
without loss of generality. 

We also require homogeneity of h — 5F/8<p, where we 
fix n c , n\,n2 and n p in the functional derivative. With 
the aid of Eq.(2.15) some calculations give 



h 
T 



1 



C' 



= tf/o(0) + ^|V0| 2 -V-CV0- 



T 



E 

=c,l,2,p 



2 

gin 



£l -E 2 



8ttT 



v V M ln(l - C), 



(2.20) 



where = dfo/d(j> and C" = dC/d4>. On the right 
hand side, the first three terms are those in the usual 
Ginzburg-Landau theory 2 ^. The last three terms arise 
from the electrostatic interaction, the solvation interac- 
tion, and the dissociation equilibration, respectively. We 
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may calculate the interface profiles and the mesophase 
profiles from the homogeneity of 

Surface tension. In the above procedure, we have min- 
imized the grand potential f2 = J druj under the charge 
neutrality J drp = 0, where the grand potential density 
is defined by 



w = f T -h<f>- Mini - p 2 n 2 , 



(2.21) 



with / T being given by Eq. (2. 7). Using Eqs. (2.16)-(2.19) 
we may eliminate [i\ and /j, 2 to obtain 

TC f 
" = h + ^r\V<t>\ 2 E 2 ~p<S> 

Z OTT 

-T J2 ni + Tv^VM^Ml-Q. (2.22) 

i=c,l,2 

Furthermore, using Eq.(2.20) we may calculate the space 
gradient of u> as 

dxk ^ dxi ( ^dxkdxe) dxk^' ^ ^ 

where d/dxu and d/dxi are the space derivatives with 
respect to the Cartesian coordinates x,y, and z. In the 
one-dimensional case, where all the quantities vary along 
the z axis, the above equation is integrated to give 



lu = TC(4>'y -p® + 



(2.24) 



where 4>' — d<j)/dz and w^, is a constant. Therefore, 
around a planar interface separating two bulk phases, 
to(z) tends to a common constant uJqo as z — > ±oo. From 
the above relation the surface tension 7 = J dz[u)(z) — 
Woo] is expressed as 



1 = 2 
T 



dz 



= I dz 



fo \ 
T~ 2- n 

i=c,l,2 

e{4>) 



VM<t> 



ln(l - C) 



C(0') 



AirT 



E A 



(2.25) 



where fo = fo — htfi — lUqo in the first line and use is made 
of J dzp<& — f dzeE 2 /4tt. In the second line the inte- 
grand consists of a positive gradient term and a negative 
electrostatic term. Similar expressions for the surface 
tension have been obtained for electrolyte and ionic 
surfactant systems^. 

Mass action law. If Eqs. (2. 12) and (2.13) are multi- 
plied, U cancels to disappear. It follows the equation of 
ionization equilibrium or the mass action equation^, 



c 



i-C 



n c = K(4>), 



(2.26) 



where K(<p) is the dissociation constant of the form, 

K (0) = Vq 1 exp[-A - (g p + g D )<f>]. (2.27) 

We may interpret A(</>) = A + (g p + g c )(f> as the 
composition-dependent dissociation energy divided by T. 



With increasing cf>, the dissociation decreases for positive 
9p + .9c and increases for negative g p + g c . If g p + g c ^> 1, 
K((j>) much decreases even for a small increase of <fi. Then 
£ and n p are related to n c as 



c 



v n p _ K 



v M (j> K + n c 



(2.28) 



These relations hold in equilibrium states, which may 
be inhomogeneous. In our theory, n p — n p (<fi,U) is a 
function of the local values of <j> and U as well as n c , n\, 
and n 2 - Here, 

^ = (1-CK, ^ = Oftnp, (2.29) 

so n p increases with increasing U at fixed <f>. 

Relations in bulk without salt. Furthermore, in a ho- 
mogeneous bulk phase with n p = n c , Eq.(2.28) yields the 
quadratic equation for n c , 



n c (n c + K) = Vq vm4>K, 
which is solved to give 

v n c 2 



C 



v M <t> VQ+T+i 

Here it is convenient to introduce 

Q(cj>) = A VM 4>/voK{4>). 



(2.30) 



(2.31) 



(2.32) 



In particular, we find £ -C 1 and n c = [vm^K/vo) 1 ^ 2 for 
Q ^> 1, while C — > 1 for Q < 1. Thus the ion density n c 
has been determined for given <fi. 

Relations in bulk with salt. As another simple situa- 
tion, we may add a salt whose cations are of the same 
species as the counterions. The cations and anions are 
both monovalent. Here the counterionsa and the cations 
from the salt are indistinguishable. Thus the sum of the 
countcrion density n c and the salt cation density n\ is 
written as n c , while the salt anion density is written as 
n 2 . The charge neutrality condition becomes n c = 



-n 2 



in the bulk. From Eq.(2.28) we obtain 

n c = Vq 1 v m <I>K/{K + n c ) + n 2 , 



(2.33) 



in the bulk phase. Then n c increases with increasing n 2 . 
We treat n 2 as an externally given constant to obtain 



= n 2 + v ~ 1 is M 4 l K/(n 2 + if), 



(2.34) 



where the second line holds in the case n 2 + K ^> 
2{v M ^K/vo) 1 ^ 2 = Q X I 2 K. With increasing n 2 , the ion- 
ized monomer density n p = n c — n 2 decreases, while n c 
increases. At high salt densities, where n 2 much exceeds 
both 2{vm4>K/ v o) 1 ^ 2 an< i K> we eventually obtain 

K<t> K . ^ 

n c = n 2 , n p = v M <§C n 2 , C — — <L (2.35) 

n 2 n 2 
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III. STRUCTURE FACTOR OF COMPOSITION 
FLUCTUATIONS 

In phase transition theories^ the order parameter fluc- 
tuations obey the equilibrium distribution oc exp(— F/T), 
where F is the Ginzburg-Landau free energy functional. 
In the present problem, the thermal flucuations of 0, n,, 
and C are assumed to obey the distribution oc exp(— F/T) 
in equilibrium, where F is the space integral of fx in 
Eq.(2.7). In the Gaussian approximation of F, we con- 
sider small plane- wave fluctuations of <fi, n,, and £ with 
wave vector q in a one-phase state. It then follows the 
mean-field expression for the structure factor S(q) of the 
composition fluctuations. It is of the form of Eq.(l.l) 
for a constant degree of ionization £ in the absence of the 
solvation interaction. Here it will be calculated including 
the solvation interaction and in the annealed case. We 
examine how it depends on the parameters g c , gi, g2, and 
g p and how it is modified by the fluctuating ionization. 



A. Gaussian approximation 

From Eq.(2.7) the fluctuation contributions to F in the 
bilinear order are written as 



SF 
T 



E [i (f+CjW + |i lwP+s ig 



E giniq4>* q + 2 (i- P C)C 2 ^ 9 



(3.1) 



where <f>q, pq, mq (i = c, l,2,p), and (q are the Fourier 
components of <j>(r), p(r), n>i(r), and £( r ): respectively. 
From Eq.(2.4) the fluctuation of the charged-monomer 
density is of the form, 



ipq 



(3.2) 



In this section <f> and denote the spatial averages, where 
Up = 53i= c i 2 Zi n i from the overall charge neutrality. 
The inhomogeneity in the dielectric constant may be ne- 
glected for small fluctuations. From the Flory-Huggins 
free energy (2.8) the coefficient f is of the form, 



v r 



vo d 2 fo 
T dd> 2 



1 



1 



1 - 



-2 X . 



(3.3) 



On the right and side of Eq.(3.1), the first term yields 
the Ornstein-Zernike structure factor of the composition 
without the coupling to the charged particles. The sec- 
ond and third terms are well-known for electrolyte sys- 
tems, leading to the Debye-Hiickel screening of the charge 
density correlation. The fourth term arises from the sol- 
vation interaction, while the fifth term from the fluctua- 
tion of ionization. 

We minimize 8F with respect to n i( j (i = c, 1, 2,p) at 
fixed 4>q to express them in the linear form ntq oc 4>q. In 



particular, pq and (q are written as 

Aq 2 



Cq = -C(l 



Ua = 



•0 



9p + 



4tv£ b A 



+ Ki 



where Uq is the Fourier component of U(r), and 



, =47r£ B [(l-C)«pH 

(i - ()g P n P 



E 

i=c,l,2 



E 

i=c,l,2 



Zigi TLi 



(3.4) 
(3.5) 

(3.6) 
(3.7) 



Here kt defined in Eq.(3.6) is a generalized Dcbye-Huckel 
wave number introduced by Raphael and Joanny^, where 
the term proportional to n p arises from dissociation and 
recombination on the polymer chains. The first relation 
(3.4) itself readily follows from the relations, 



dp 



-eA, 



dp 
dU 



'4arE 



(3.8) 



where p(r) is regarded as a function of the local values of 
U(r) and <j>(r) as discussed around Eqs.(2.28) and (2.29). 
From Eq.(3.5) we find (q —* as £ — > or 1. We can see 
that the last term on the right hand side of Eq.(3.1) is 
negligible as £ — ► 1- 

After elimination of rnq (i = c, 1, 2,p) the free energy 
change SF is expressed as SF = ^2qT\(j>q\ 2 / S(q), where 
S(q) is the composition structure factor calculated as 

1 „ o 4:Tr£ B A 2 



S(q) 



■Of 



f + r s + CYp^T + C T 



1 



lp K T 



(3.9) 



where r s is the shift of r arising from the solvation inter- 
action given by 



= pr 1 



9p 



a 



Qgl\ n P 



E 

i=c,l,2 



(3.10) 



As <f> — » we have r s = 2g p f. For large ^> 1 the 
other negative terms (oc g\) can be significant, leading to 
r s < 0. In the second line of Eq.(3.9) 7 P is defined by 



7 P 



C Krp 



(Att£ b C)-^ 2 A 



Op + Ei= Cl i,2 



(3.11) 



where (^kI b C)- x I 2 = [3a0(l - 4>)/tt£ b ] 1/2 if use is made 
of Eq.(2.9). Note that r s and consist of the terms 
proportional to the charge densities, while "f p depends 
on their ratios. In particular, if n c = n p and n\ = n<i 
in the monovalent case, 7 P depends on the density ratio 
R = ni/n p as 



7p 



(4^ B C) 



-1/2 



2 - C + 2R 



T - (1- C)9p+9c + (.91 



92)R ■ 

(3.12) 
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Here R — for salt-free polyelectrolytes and R is in- 
creased with increasing the salt density. As R — > oo the 
above formula tends to that for neutral polymer solutions 
(low-molecular- weight binary mixtures for N = 1) with 
salt. We note that 7 P can even be negative depending 
on the solvation terms in the brackets on the right hand 
side of Eq. (3. 12). 

B. Macrophase and microphase separation 

From the second line of Eq.(3.9) we obtain the small-q 
expansion S(q)^ 1 = ao + o-iq 2 + 02<? 4 ■ ■ • with 

at = C(l - 7p 2 ), (3.13) 

so we encounter the following two cases, (i) If |7 P | < 1, 
S{q) is maximum at q = and we predict the usual phase 
transition with increasing \ m the mean field theory. The 
spinodal \ = Xsp is given by 

2x sp ~ + T ^+r s + C 7p '4- (3.14) 

Macroscopic phase separation occurs for x > Xsp- On 
the right hand side of Eq.(3.14), r s = and the last term 
becomes n p /(l — Qcp 2 > for g t = 0, but the sum of 
the last two terms can be much altered for \gi\(j> > 1. 
For example, without salt and in the monovalent case, 
it is equal to [1/0 2 + (g 2 + 2g c /4>]n p /(l - () - 2g p g c n p . 
(ii) If |7 P | > 1, S(q) has a peak at an intermediate wave 
number q m given by 

ql = {^l B /C)^\A\-^ 

= 4(|7pI"1). (3-15) 

Thus q m — > as kt — > or as | 7p | — > 1. The spinodal is 
given by S^m) -1 = or by 

2Xs P = + + r s + C4(2\ lp \ - 1). (3.16) 

Microphase separation should be triggered for \ > Xsp- 
Since 7 P depends only on the ion density ratio R as in 
Eq.(3.12), the criterion 7 p > 1 remains unchanged how- 
ever small the densities of the charged particles are. Of 
course, mesophase formation is well-defined only when 
q m in Eq.(3.15) is much larger than the inverse of the 
system length. In the following we discuss some special 
cases. 

Polyelectrolytes without solvation interaction. For gi = 
0, we have r s = and A — n p /<p = v^ 1 / as in the 
previous theories of weakly ionized polyelectrolyte a 5 ' 9 ' 10 . 
In the monovalent case Eq.(3.12) gives 

1 P = l/[(4 7 r^C) 1 / 2 (2-C + 2i?)0] 

= [3a^AMW)] 1/2 /(2-C + 2i?), (3.17) 

where R — m/n p and use is made of Eq.(2.9) in the sec- 
ond line. Here 7p is largest without salt (R = 0) and 



decreases with increasing R. In accord with this pre- 
dicted salt effect, Braun et aL 1 - observed a mesophase at 
low salt concentration and a macroscopic phase separa- 
tion at high salt concentration. However, for large |^|, 
the solvation interaction comes into play. Thus more ex- 
periments are desirable to detect it in polyelectrolytes 
with water- like solvent. 

Binary mixtures with solvation interaction. Without 
ionization or for / = 0, we have n p = n c = 0. Here we 
describe a neutral binary mixture with salt. Further for 
a monovalent salt, we have A = (g\ — g2)n\ and 

1 P = (gi-g2)/{lM B C) 1 f\ (3.18) 
r s = -(gl+g^ni. (3.19) 

From neutron scattering, Sadakane et al^ found peri- 
odic structures in a near-critical, low-molecular-weight 
mixture of D2O and trimethylpyridine (3MP) containing 
sodium tetrarphenylborate (NaBPlu). Their salt is com- 
posed of strongly hydrophilic cation Na + and strongly 
hydrophobic anion BPli4~. If their data are interpreted 
in our theoretical scheme, we expect g\ 3> 1, 32 <C — 1, 
and 7p > 1 for their system. Such ion pairs are antag- 
onistic to each other, undergoing microphase separation 
where the solvent composition is inhomogeneous. It is 
worth noting that the scattering amplitude was maxi- 
mum at q = upon addition of NaCl, KC1, etc. in the 
same mixture D2O-3MP 35 . These salts consist of hy- 
drophilic anions and cations with g\ *~ g%, so 7 P should 
be much smaller leading to no mesophase formation. It 
is also striking that the coexistence curve is much shifted 
as AT (x r s (x gf for large |<jfj| with increasing the salt 
density even for hydrophilic ion pairs, which is consistent 
with a number of previous experiment a 36 ' 37 ! 38 . 

Neutral polymer solutions with salt. For very weak ion- 
ization but in the presence of salt, 7 P is expressed as in 
Eq.(3.18). In our theory, even polymer solutions consist- 
ing of neutral polymers and a polar solvent can exhibit 
a charge-density wave phase for \g± — 52I S> 1< Hakim et 
a h 21 ! 22 found a broad peak at an intermediate wave num- 
ber in the scattering amplitude in (neutral) polyethylene- 
oxide (PEO) solutions with methanol and with acetoni- 
trile by adding a small amount of salt KI. They ascribed 
the origin of the peak to binding of potassium ions to 
PEO chains. In our theory such a peak can arise for a 
sufficiently large g\ — gi- Remarkably, the peak disap- 
peared if the solvent was water, which indicates sensitive 
dependence of the molecular interactions on the solvent 
species. Thus more experiments should be informative, 
where use of antagonistic ion pairs will yield more drastic 
effects, leading to mesophases. 



IV. ONE-DIMENSIONAL PROFILES AND 
NUMERICAL RESULTS 

We give numerical results of one-dimensional profiles 
in equilibrium. All the quantities are assumed to vary 
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along the z axis. For each given 4>(z), we need to solve 
the nonlinear Poisson equation (2.10) numerically. The 
charge density is expressed as 



(a) A =5, 9 = 1, 9 r =4 



E 

i=c,l,2 



Z i n l (z) 



e.v M <t>(z) 



l + n c {z)/K(z) 



(4.1) 



where rii(z) oc exp[—gi(f>(z) — ZiU(z)] from Eqs.(2.17) and 
(2.19) and K(z) depends on <f>(z) as in Eq.(2.27). We 
then seek <f>(z) self-consistently from the homogeneity of 
the right hand side of Eq.(2.20). 

As in electrolytes, there appears a difference of the elec- 
tric potential A<I> across an interface, which is called the 
Galvani potential difference^. Because there are many 
parameters, we will set % = 1, iV = 20, E\ — — 0.9eo, an d 
Ib = e 2 /sqT = 8a/ir in all the following examples. Then 
the dielectric constant of the solvent is 10 times larger 
than that of the polymer. The space will be measured in 
units of the molecular size a 



1/3 



A. Interface without salt 

First we suppose coexistence of two salt-free phases, 
separated by a planar interface. Hereafter the quantities 
with the subscript a (the subscript 0) denote the bulk 
values in the polymer-rich (solvent-rich) phase attained 
as z — > — oo (as z — > oo). Namely, we write n ca = 
n c (— oo) and n c p — n c (oo). Here we obtain n ca or n c p 
from Eq.(2.29) and 



v Co 



n c (3 



exp(-g c A0 - AU) 



(4.2) 



from Eq.(2.17). Hereafter A<j> 
tion difference and 



<j>p is the composi- 




40 45 50 55 60 35 40 45 50 55 60 

z/a z/a 



(b) A =6, 9 p =2, g c =-6 




25 30 35 40 45 50 25 30 35 40 45 50 

z/a z/a 

FIG. 1: Interface profiles in the salt-free case for (a) Ao = 5, 
g p = 1, and g c = 4 (top) and for (b) Ao = 8, g p = 2, and 
<7 C = —6 (bottom). Polymer volume fraction <j>(z), normalized 
potential <&(z)/<&m, and degree of ionization £(z) (left), and 
normalized charge densities von c (z) and von p (z) (right). The 
other parameters are common as \ = 1, iV = 20, vm = 0.1, 
ei = — 0.9eo, and Ib = 8a/-K. Here $(z) is measured from its 
minimum, and $Af(= 2.67T/e in (a) and 6.86T/e in (b)) is 
the difference of its maximum and minimum. 



AU = eA$/T = e($« - ®p)/T 



(4.3) 



is the normalized Galvani potential difference. In terms 
of Q a = Av M (j) a /v K a and Qp = Av M 4>{3/v Kp, we ob- 
tain 



AU = g p A(f> + ln 







i - 1 



y/Q a + 1-1 

My/Qp + l + l)' 



(4.4) 
(4.5) 



In particular, if Q a ^> 1 and > 1 (or ( Q < 1 and 
(fj < 1), we obtain AU = (g p - g c )Acf>/2 + \n{^>p/<j) a ) 
and n cl3 /n ca = exp[(g p + 5 c )A0/2](0 (3 /0 Q ) 1 / 2 . 

The interface profiles are extremely varied, sensitively 
depending on the molecular interaction parameters, Ao, 
g p , and g c . As such examples, in Fig.l, we show salt-free 
interface profiles for (a) Ao = 5, g p = 1, and g c = 4 and 
(b) Ao = 8, g p = 2, and g c = —6. In the a and /3 regions, 
the degree of ionization ( is 0.071 and 0.65 in (a) and is 
0.24 and 0.51 in (b), respectively. The charge densities 



are multiplied by v = a 3 . The AU = e($ Q - ®p)/T is 
-2.67 in (a) and 0.099 in (b). Interestingly, in (b), U(z) 
exhibits a deep minimum with C/ m i n = Up — 5.98 at the 
interface position. We can see appearance of the charge 
density n c — n p around the interface, resulting in an elec- 
tric double layer. The counterion density is shifted to 
the (3 region in (a) because of positive g c and to the a 
region in (b) because of negative g c . The parameter j p 
in Eq.(3.11) is 0.75 in (a) and 0.20 in (b) in the a region, 
ensuring the stability of the a region. From Eq.(2.25) the 
surface tension 7 is calculated as 7 = 0.0175T/a 2 in (a) 
and as 0.0556T/a 2 in (b), while we obtain 7 = 0.050T/a 2 
without ions at the same x = 1. In (a) the surface ten- 
sion 7 is largely decreased because the electrostatic term 
in Eq.(2.25) is increased due to the formation of a large 
electric double layer. In (b), on the contrary, it is in- 
creased by 10% due to depletion of the charged particles 
from the interface^. 
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(a) A =5, g - 1, 9 C =4, g 2 =2 
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(b) a =5, o - 1, g c =4, g 2 
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FIG. 2: Interface profiles with mobile cations and anions 
for (a) A = 5, g p = 1, g c = 4, and g 2 = 2 (top), for (b) 
Ao = 5, g p = 1, g c = 4, and 32 = —8 (middle), and for (c) 
Ao = 8, (j p = 2, g c = —6, and g 2 = 2 (bottom). Curves repre- 
sent </>(z), &(z)/$m, and £(2) (left) and uoc c (z), voCp(z), and 
«oC2(z) (right) with won-2/3 = 0.002. Here $(z) is measured 
from its minimum, and $m is the difference of its maximum 
and minimum. The other parameters are the same as those 
in Fig.l. 



B. Interface with salt 

With addition of salt, interface profiles are even more 
complex. For simplicity, we consider a salt whose cations 
are of the same species as the counterions from the poly- 
mer. This is the example discussed at the end of Sec- 
tion 2. The free energy density in this case is still given 



by Eq.(2.7) if we set n\ = there. The densities of 
the mobile cations, the mobile anions, and the charged 
monomers are written as n c (z), 712(2), and n p (z), respec- 
tively. All the charged particles are monovalent. We 
treat 712/3 in the solvent-rich (3 region as a control param- 
eter, which is the salt density added in the (3 region. The 
densities 772(2) and n c (z) depend on z as 



772(2) = 77 2 /3 exp[g 2 ((j)p - H z )) 

n c (z) = 77 c/3 exp[g c (0/3 - <p{z)) - 



Up + U(z)}, 
U -U(z)}. (4.6) 



We write S = 772/3 / Kp to avoid the cumbersome notation 
in the following. Use of the first line of Eq.(2.34) gives 



77 c /3 



2 



S -1 + ^(S + l) 2 + Q/3 



(4.7) 



We use Eq.(2.33) in the a region to determine n ca . Since 
n ca n 2a /K a = 77 c/3 n 2 p/Kp from Eq.(4.5), n ca obeys 

(t7 cq - Sn c pK a /n ca )(l + n ca /K a ) = v^VM&a- (4.8) 

In particular, if Ti CQ Kq, ot if C,a. < ^ I5 W6 find 

77 CQ S Ay 2 (577 c/3 + v^v M 4> a ) 1/2 , (4.9) 



AU S g p A(f> + - In 



(4.10) 



Sn c p/K a + Q a 

In Fig. 2, we show profiles of 4>(z), Q(z), C( z ), n c (z), 
772(2), and n p {z) with 772,9 = 0.002tj^~ 1 held fixed for three 
cases: (a) A = 5, g p = 1, g c = 4, and g 2 = 2; (b) 
A = 5, g p = l,g c = 4, and g 2 = -8; (c) A = 8, g p = 2, 
g c = —6, and g 2 — 2. The normalized potential difference 
AU is -1.61 in (a), -2.54 in (b), and 1.88 in (c), while 
the normalized surface tension a 2 7/T is given by 0.0358 
in (a), 0.0233 in (b), and 0.0554 in (c). In (a) the cations 
and anions are both repelled from the interface into the 
region because of positive g c and g 2 . In (b) , owing to g 2 = 
—8, the anions become richer in the a region than in the 
[3 region, since they are strongly attracted to the polymer 
chains. There, despite g c — 4, the cations become also 
richer in the a region to satisfy the charge neutrality. In 
(c), we set g c — —6, so the cations are strongly attracted 
to the polymer chains. The cations are also attracted 
into the a region. The electric dipole due to the double 
layer is reversed, yielding AU = U a — Up > 0. 



C. Periodic state 

With varying the temperature (or \) 1 the average com- 
position (0), the amount of salt, there can emerge a 
number of mesophases sensitively depending on the var- 
ious molecular parameters (<?i, Ao, and Vm)- In Fig. 3, 
we show an example of a one-dimensional periodic state 
without salt. Here vm is set equal to 0.5 and the charge 
densities are much increased. In this case, the degree 
of segregation and the charge heterogeneities are much 
milder than in the cases in Fig.l. Since the charge density 
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FIG. 3: Periodic profiles in a salt-free mesophase with vm = 
0.5, A = 5, g p — 1, and g c = 4. Top: 4>{z). <&(z)/<&m with 
<E>ju = 0.86T/e, and ((z). Bottom: von c (z) and uon p (z). 



p(z) everywhere remains small, Eq.(2.31) locally holds 
approximately. Thus, 



n c (z) = 



VM<j>(z) 



v/Q(z) + l/4+l/2 



(4.11) 



where Q(z) is defined by Eq.(2.32) and n p (z) is deter- 
mined by Eq.(2.26). The deviation of n c (z) from the 
right hand side of Eq.(4.10) gives $(z) as a first approxi- 
mation. The above relation approximately holds slightly 
below the transition from a disordered state to a charge- 
density- wave state. 



V. SUMMARY AND REMARKS 

The charge distributions in polyelectrolytes are ex- 
tremely complex around interfaces and in mesoscopic 
states, sensitively depending on the molecular (solvation) 
interaction and the dissociation process. Our contin- 
uum theory takes account of these effects in the simplest 
manner, though it should be inaccurate on the angstrom 
scale. 

Our main results are as follows. In Section 2, we have 
presented a continuum theory accounting for the sol- 
vation effect and the dissociation equilibrium. The de- 
gree of ionization is rather a dynamic variable dependent 



on the local electric potential. The surface tension in 
Eq.(2.25) consists of a positive gradient term and a neg- 
ative electrostatic term, where the latter is significant as 
the electric double layer is enhanced as in the top pan- 
els in Figs.l and 2. In Section 3, we have calculated the 
structure factor S(q) of the composition fluctuations for 
polymer solutions with the solvation and ionization ef- 
fects included. When the parameter 7 P in Eq.(3.11) or 
(3.12) exceeds unity in a one-phase state, a peak can ap- 
pear at an intermediate wave number in S(q). In such 
systems, a mesoscopic phase can emerge as the interac- 
tion parameter \ is increased. In Section 4, we have nu- 
merically solved the equilibrium equations in Section 2 to 
examine one-dimensional profiles of interfaces with and 
without salt and in a periodic state without salt. Though 
very preliminary and fragmentary, our numerical results 
demonstrate dramatic influence of the solvation interac- 
tion on inhomogeneous structures. 

A number of complex effects should be further taken 
into account to properly describe polyelectrolytes. (i) 
Under Eq.(2.1) we have neglected the effect of the elec- 
trostatic interaction on the chain conformations to use 
the Flory-Huggins free energy density. Hence, as argued 
in the literature 2 -! 3 -, our theory cannot be justified at small 
4>. (ii) We have neglected the image interaction, which 
is known to increase the surface tension of a water-air 
interface at low ion densities^. Generally, it arises from 
a deformation of the self-interaction of ions due to inho- 
mogeneous £~. So it originates from the discrete nature 
of ions and is not accounted for in the free energy (2.7), 
where the electric field $ is produced by the smoothly 
coarse-grained charge density p. (iii) At sufficiently low 
ion densities, we have neglected the ionic correlations giv- 
ing rise to ion dipoles and cluster o 12 i 28 and a free en- 
ergy contribution due to the charge density fluctuations 
(oc k 3 ) 25 ' 26 . (iv) The phase diagram of mesophases in- 
cluding the solvation interaction remains unknown. To 
describe the mesophases we need to perform nonlinear 
analysis and computer simulations, (v) We have assumed 
a one-component water-like solvent. For two-component 
solvents we may predict a variety of effects. For example, 
with addition a small amount of water in a less polar sol- 
vent, hydration shells are formed around the hydrophilic 
ionized monomers and the counterion a 8 ' 33 , leading to an 
increase of the degree of ionization. For not small water 
concentrations, we may well predict formation of com- 
plex mesoscopic structures mediated by the Coulombic 
and solvation interactions. 
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